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ABSTRACT: Two-temperature advection dominated accretion flow (ADAF) or hot ion 
tori (HIT) models help explain low luminosity stellar and galactic accreting sources and 
may complement observational support for black holes in nature. But low radiative effi- 
ciencies demand that ions receive a fraction r\ > 99% of energy dissipated in the turbulent 
accretion. The r\ depends on the ratio of particle to magnetic pressure. If compressive 
modes of dissipation, like magnetic mirroring, dominate incompressible modes, then even 
when the pressure ratio is 0(1), the required large rj can be attained. However the relative 
importance of compressive vs. incompressible modes is hard to estimate. The larger up 
in the turbulent cascade, the more compressible the turbulence. The relevant length scale 
for particle energization can be determined by equating the dominant eddy turnover time 
to the time for which an energy equal to that in the turbulence can be drained. Based 
on the large scales resulting from this estimate, it is suggested that compressive mirroring 
may be important. Also, regardless of the precise rj or dissipation mechanism, non-thermal 
protons seem natural in two- temperature discs because all dissipation mechanisms, and 
the use of an isotopic pressure, require wave-particle resonances that operate only on a 
subset of the particles. Finally, it is briefly mentioned how mirroring may help to generate 
an ADAF or HIT in the first place. 
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1. Introduction 



Magnetized accretion discs have become the most convincing physical paradigm to 
explain emission from the central engines of active galactic nuclei (AGN) and X-ray binary 
sources (Frank et al. 1992). The observed radiation comes from the energy dissipation 
required to maintain steady accretion of material onto the central object. As molecu- 
lar viscosity is incapable of providing the required accretion rates, turbulent viscosity is 
necessary. For thin discs, this can be generated by shear and magnetic fields (Balbus & 
Hawley 1991). For thick discs, something similar may occur, though in this case angular 
momentum transport may ultimately require a global approach. 

Nevertheless, in complement to thin disc solutions for sources requiring high radiation 
efficiency accretion (e.g. Frank et al. 1991), two-temperature thick advection dominated 
flows (ADAFs) or hot ion tori (HIT) (e.g. Ichimaru 1977; Paczynski & Bitsnovatyi-Kogan 
1981; Rees et al. 1982; Narayan & Yi 1995) have received much attention in an effort to 
explain sources requiring a low radiation efficiency. Here the ions are assumed to receive 
the energy dissipated by the steady accretion without having enough time to transfer their 
energy to the cooler electrons before falling onto the central object. Some or most of the 
dissipated energy is advected, not radiated, as it would have been if electrons received all 
of the dissipated energy. Such models have been at least partially successful in explaining 
quiescent galactic centres (Rees 1982; Narayan et al. 1995; Mahadevan 1998; Fabian et al. 
1995, but see DiMatteo et al. 1998) and stellar X-ray binary systems (Narayan et al. 1996) 
with radiative efficiencies < 1/100 that of thin disc solutions. When the central object is 
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a black hole, the advected energy is lost forever rather than re-radiated as it would be for 
a neutron star. Precisely such observed differences between corresponding X-ray binary 
systems have been purported to provide evidence for black hole horizons (Narayan et. al 
1997). 

There have been only a handful of other papers addressing how the viscous dissipa- 
tion might energize particles in accretion flows (Gruzinov 1997; Quataert 1997; Quataert 
& Gruzinov 1998) and none addressing thermal vs. non-thermal particle distributions. 
Both of these issues are extremely important for ADAF type models because: 1) a two- 
temperature solution is not sufficient to explain a low radiation efficiency and 2) inter- 
pretation of observations of the Galactic centre suggest that the protons are non-thermal 
when an ADAF model is employed (Mahadevan 1998). The potential catastrophe for 
ADAF/HIT models, if electrons are preferentially energized over protons, was partially 
explored in (Bitsnovatyi-Kogan & Lovelace 1997). Because electrons cool much faster than 
ions, even if 1/2 of the dissipated energy went into electrons, a two-temperature solution 
would still result. But in this case, 50% of the dissipated energy would be radiated — far 
too much to explain low luminosity sources. More explicitly, I define qt as the magnitude 
of the energy density input rate into particles. Then q t = g+ + g+ = rjq t + (1 — r})q t 
where q p and q+ are the magnitudes of energy input rate into protons and electrons, 
and r\ is the fraction of q t that goes into protons. In the steady state, energy loss rates 
are equal to energy gain rates so that when advection is included, we have for the pro- 
tons q+ = r\q t = q~ - q pe = fr\q t + (1 - f)rjq t , where q a = fr\q t is the rate associated 
with advection, q pe is the rate of transfer from ions to electrons, and / is the fraction of 
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the proton energy loss rate associated with advection. For the electrons, we thus have 
q+ = (l — r])q t + q pe = Qt — Qa = <7t(l — fv)- Since q+ = q~ ', where q~ is the luminosity den- 
sity, the quantity (1 — frj) must be << 1 to explain quiescent sources. Standard treatments 
(e.g. Rees et al. 1982; Narayan & Yi 1995) assume r\ = 1 so that 1% radiative efficiency 
would correspond to / = 0.99. The important questions are: (A) When can rj > 0.99 be 
justified? (B) Are protons non-thermal? and (C) Is there a faster than Coulomb coupling 
between electrons and protons (e.g. Begelman & Chiueh 1988) that destroys the ADAF 
solution? I will address (A) and (B) here. 

Proton energization by incompressible (Quataert 1997; Gruzinov 1997; Quataert & 
Gruzinov 1998) and/or compressive modes of dissipation (discussed herein) both depend 
on the ratio of particle to magnetic pressure. Section 2 addresses the relation between 
magnetic, turbulent, and particle energy densities in ADAFs, relating them to the viscosity 
parameter a. Section 3 discusses the threat of electron runaway. Section 4 employs a very 
physical approach to acceleration by magnetic mirroring and derives the time scale for 
particle energy doubling for two distinct limits of the average particle speed. Section 5 
discusses energization in ADAFs, first addressing why protons are likely to be non-thermal, 
regardless of the acceleration mode. The mirroring results are then specifically applied to 
ADAFs and the scale in the turbulent cascade where mirroring is favored is estimated. 
Because larger scales are significantly compressive, and the resulting derived scale can be 
large, mirroring may be important. Magnetic mirroring type processes can favor protons to 
the extent required for a wider range of average particle to magnetic pressure ratios than 
found by Quataert & Gruzinov (1998) from dissipation of incompressible Alfven waves, 
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but the relative fractions of incompressible vs. compressible modes of dissipation are hard 
to determine. The possibility that mirroring may help provide a thermal instability which 
initially forms an ADAF is briefly addressed. 

2. Relation Between Viscosity Parameter and Pressure Ratio 

The standard parameterization of accretion disc turbulent viscosity for thin discs is 
(Shakura & Sunayev 1973) 

v T = aC s H ~ Vt/t/3, (1) 

where a is the viscosity parameter and C s and H are the sound speed and disc height, 
and Vt and It are the outer (i.e. dominant energy containing) random (turbulent) flow 
speed and scale. For thin discs, the scale of the turbulence is always much less than 
the radius of the disc. Turbulent kinetic and magnetic energy densities rapidly approach 
equipartition from field line stretching. A magneto-rotational shearing instability (Balbus 
& Hawley 1991) likely drives 3-D MHD turbulence. In a steady state, dissipation of the 
turbulence into particles combats the symbiotic growth of magnetic and kinetic turbulent 
energy. Both magnetic and kinetic turbulent energies incur a decaying power-law energy 
spectrum (like Kolmogorov (1941) or Kraichnan (1965)) with the largest scales of the 
turbulence containing the most energy. Since the sound speed is constant on all scales, 
the largest scales are thus the most compressive. For thin discs, the outer turbulent scale 
is significantly smaller than the disc radius, but for standard ADAFs, there is not such a 
strong scale separation (Blackman 1998). 

For thin discs, we can derive a 1-to-l link between a and = Vj[/Cg = 6(1 — /3 a ), 
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where Va is the Alfven speed, and (3 a is used in ADAF modeling. In the steady-state, 
the largest eddy turnover time tr = It/Vt must equal the shearing instability growth 
time driving the turbulence, that is: tr — R/Vcf,, where the rotation speed ~ Vk, the 
Keplerian speed. Since Va ~ Vr in the saturated state from field line stretching (e.g. 
Parker 1979; Blackman 1998) and C s /Vk = H/R from hydrostatic equilibrium, we have 
v = aC s H ~ V T l T /3 ~ VKR/ZVt) = V^(H/3C S ) which implies that 

a ~ 2(1 - p a )(V K /V4,) ~ 2(1 - a ) (2) 

for thin discs. This result basically agrees with numerical simulation results (e.g. Stone 
et. al 1996). 

For ADAFs, the ratio of Vk/V^ can be so low that (2) is inappropriate: in this case 
the resulting eddy scale implied by the relation would be larger than H ~ R. We can 
instead obtain an upper limit on a for ADAFs that comes from the constraint 

l T <H ~R. (3) 

Then from (1), the definition of [3 a , and Vr = Va from field line stretching, we find 

a<(2/3 1 /2)(i_ /3a )V2. (4) 

showing that j3 a and a are not independent. 

3. On Electron Runaway 

Bistnoyati-Kogan & Lovelace (1997) pose an interesting question: why can't direct 
acceleration from electric fields drain energy into electrons, destroying ADAFs? I address 
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this here. First, note the generalized Ohm's Law (e.g. Scudder 1986) E = — V e /c x B + 
cr _1 J — m e (V e - VV e )/e — VP e /(en e ) where B is the magnetic field, J is the current density, 
P e is the electron pressure, V e is the bulk electron velocity, a -1 is the resistivity, m e is the 
electron mass, n e is the density, and e is the electric charge. For the plasmas of interest, 
a characteristic magnitude of E parallel to B is given by the second last term, and using 
V e <V T ~V A < C s , we have ~ k b T e /(5l e)~2x 10- 14 (T e /10 9 K)(<W/10 13 cm)-\ 
where 51 is the gradient length, T e is the electron temperature and k b is the Boltzmann 
constant. 

For £j| to produce ER, it would have to exceed the Dreicer electric field (Dre- 
icer 1962; Holman 1985; Bitsnovatyi-Kogan & Lovelace 1997) E D = e 2 lnT/X 2 d = 1.8 x 
10- 7 (/nr/20)(n e /10 12 cm- 3 ) 1 / 2 (Te/K^K^St- Volt/cm, using Ad ~ 6.65(T e 1/2 /n e )cm, 
for the Debye length. Whether E\\ > E D depends on the size of 51. For AGN, E D 
is only exceeded on scales 10 6 times smaller than the turbulent outer scale. For stel- 
lar size X-ray binary ADAF systems, the outer scale is ~ 10 7 cm, so in principle ER is 
possible throughout the flow. But the accelerated electrons can never produce a cur- 
rent which induces a magnetic field in excess of the inferred ambient field. This gives 
an upper limit (Holman 1985) to the size of field gradients that generate ER, namely 
51 < 8(S/10 4 G)(n e /10 12 cm- 3 )- 1 (T e /10 9 K)- 1 / 2 ( J E; D / J E || )cm. For all relevant accretion 
discs, this scale in the cascade is always way below that at which magnetic mirroring, em- 
ployed in the next section, could have already drained most of the energy in the cascade. 
Nevertheless, if mirroring is not important, or equivalently, if a significant component of 
the turbulence cascades to incompressible scales before draining, then the cascade may 
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proceed down to this scale where ER, or other (e.g. Quataert & Gruzinov 1998) electron 
energization processes may be important. 

However, the more extreme ER of Bistnovati-Kogan & Lovelace (1997) is not likely. 
They employ the mean electric field, obtained by coarsely averaging E = (E) + Ex over 
the turbulent scales It- The dominant terms in this mean Ohm's law are then (E) ~ 
— (V T x B T ) — (V) x (B) where the turbulent EMF (Parker 1979) is, in kinematic theory, 
(Vt x By) = (ctd/c)(B) — /3dV x (B) where ad is a pseudo-scalar helicity, and (3d is the 
turbulent diffusivity. A representive magnitude of (E)\\ using |a d | ~ Vr/3, is (£j|) ~ 
V T (B)/3c. Assuming that (B) ~ B, then (E\\) ~ 3 x 10 3 (y T /10 10 cm/sec)(S/10 4 G)St- 
Volt/cm. Thus (E) » E D and one might be tempted to conclude that runaway electron 
acceleration is extreme. But particles do not actually see (E), since the average is only 
defined on scales larger than l T . Thus extreme ER should not occur. 

4. Energization by Magnetic Mirroring 

4.1 Basic Physical Picture 

Fermi energization, or magnetic mirroring of particles off magnetic compressions 
(Fermi 1949; Spitzer 1962) provides a means of dissipation of compressive turbulent energy 
into particles. Condsider a field, Bt, which represents the field on the largest turbulent 
scale, superimposed on which is a smaller scale compression, so the total field in the com- 
pression is B = Bt + SB. These compressions travel along field lines at speeds ~ Va, 
and can transfer energy to particles: consider what happens as the particle traveling along 
Bt interacts with the compression Bt + SB. Since the magnetic force is perpendicular 
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to the particle velocity, as long as the magnetic gradient scale >> particle gyro-radius 
(adiabatic approximation), the particle's angular momentum and energy are conserved in 
the frame of a magnetic compression at rest. Denoting quantities in this frame by a prime 
and working in the non-relativistic limit, the energy and angular momentum magnitudes 
are given by u' = mv' 2 /2 and j' = mv' ± r g = m 2 cv±_' 2 /(eB), where m is the particle mass, 
v', v' ± are the total speed and speed perpendicular to the field and r g is the gyro-radius. 
The constancy of both u' and j' implies that v' 2 /B = v' 2 sin 2 '(f)' j B is also constant. Thus, 
because v' is constant, sin 2 (p' oc B = B T + 5B : or 

si n 2 (j)' = sin 2 cj)' T (B T + 5B)/B T . (5) 

When sine/)' = 1, the particle reflects. Thus there exists a minimum pitch angle the 
particle must have with respect to the ambient Bt such that it can reflect upon entering 
the compression. This is given by 

sin 2( pT,imn = sin 2 (f)' min = B T /(B T + SB) (6) 

(see also figure 1). In the lab frame, the moving compression then boosts the velocity 
component of a given particle parallel to Bt- For energy to be gained from repeated 
reflections, the boost must be rapidly isotropized by particle generated waves (Eilek & 
Hughes 1991; Larosa et al. 1996) as discussed further in section 5. Assuming isotropy in 
the lab frame, in the frame of a moving magnetic fluctuation the velocity distribution is 
centered around ~ Va- The minimum angle for mirroring by a magnetic compression then 
gives a minimum speed that particles need to reflect: 

Vmin = V A sin4>' min = V A B T /2 /(B T + 5B) 1 / 2 ~ V A for 5B/B T < 1, (7) 
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as seen in figure 1. 

4.2 Time Scale for Energy Doubling 

Different regions will have Bt aligned in different directions but consider Bt in one 
region of size It in which Bt is assumed constant. Following (Larosa et al. 1996), define 
r r = U r (dU r /dt)~ 1 = N5t as the time scale for the average reflected particle energy U r to 
double, where N is the number of required reflections and St is the time between reflections. 
Since only a fraction F are reflected, the energization time averaged for all particles is then 

r = U/dU/dt = (5t/ F)(NU/U r ), (8) 

where U is the average particle energy averaged over all particles, we need to estimate 
N, St, F, U/U r . 

To understand the role of F, two limits not usually distinguished must be considered 
separately: For particles with speeds v m i n < v < Va, all reflections are "head-on" since 
the particles can never catch up to the fluctuations which move at speeds along the field 
lines ~ Va (see figure 1). For v » Va, there are both "catch-up" and head-on reflections. 
Let LI and L2 label separate regimes where the average particle speed, v ave , satisfies for 
LI: v ave » Va ~ v min and L2: v ave ~ v min ~ Va- That v min ~ Va follows from the 
assumption that SB < Bt in (7). LI is appropriate for electrons in a thermal equilibrium 
system whose magnetic pressure is not dominant. L2 is appropriate for protons in a 
thermal equilibrium system, or for both protons and electrons in plasma of j3 p ~ 1 with 
the ratio of proton to electron temperature T p /T e > 1000 — like ADAFs or HIT. Define 
the corresponding energization times tli and tl2 for the two limits. Limit LI is the 
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standard "stochastic Fermi acceleration" limit, and the energization time for LI derived 
below simply, agrees with other treatments (e.g. Miller 1985, Melrose 1986). The limit L2 
produces a different formula. 

To proceed further, I compute velocity moments of reflected particles in the two cases 
LI and L2 for thermal and non-thermal distributions. This is necessary for computing F. 
For a power-law distribution, dg nt = (A — l)(v/vo)~ x d(v/vo), where dg nt /d(v /vq) is the 
distribution function, vq is the lower cutoff on the power-law, and A > 3 will be assumed 
(to avoid the appearance of logarithms) . Integrating over the reflected particles gives 



f 



dg n t = (Vmin/Vo) 1 \ (9) 



if vq < Vjnin. If vq > v m i n , then vq replaces v m i„ as the lower integral bound. The 
average velocity of reflected particles is then v r ^ ave \ nt = (A — l)vo(v m i n /v ) 2 ~ x /(X — 2), 
when v < v min . In the thermal case, dg t h = (4:/n 1/2 )(v 2 /vl ve )Exp[-v 2 /vl ve ]d(v/v ave ), 
so that 

/oo 
dg th = 1 - Erf[v min /v ave ] + (v min /v ave )Exp[-Vm in /vl ve ]. (10) 
- min 

I now use (9) and (10) to determine F, the fraction of particles that reflect. Generally, 
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fdg, (11) 



III 171 



where / is the fraction of particles that reflect at a given speed. The / is the "area" of 
the sphere (see Figure 1) corresponding to particles which can be reflected, divided by the 
total area, Anv 2 : 

v±[l + (dv ± /dv\\) 2 ] 1/2 dv\\/47rv 2 = 27r^[^ N ]^[+ /4ttv 2 , (12) 

Hi- 
ll 



where ||(_L) indicates parallel (perpendicular) to B Tl and the second equality comes from 
using the equation for the circle centered at Va for the particle speed, v\ = v 2 — (v\ \ — Va) 2 - 
The bounds U|| are determined by finding the abscissa values at which the line defining 
tfimin intersects (Figure 1) the circle defined by v ave . Setting the equation for the lines, 
v ± = v f\tanS >J 4>min, equal to that of the circle gives 

v\\ ± = V A coscj)' min ± cos(j)' min (v 2 ave - VAsin 2 (j)' min ) 1/2 . (13) 

For LI (i.e. v ave » V A ), using (13) in (12) gives / ~ cos<\>' min ~ {5B/B T ) 1 ' 2 ~ 
2(2Tcv 2 cos(f)' miri )/4nv 2 = coscf)' . This can be pulled out of the integral in (11). Then 
because v m i n /v ave « 1 in this limit, 

F L1 = f = cos(P' min , (14) 

for both the non-thermal and thermal cases. 

For L2, (v ave ~ Vmin) (12) and (13) give / <cos(j)' min . This can be pulled out of the 
integral in (11). For the non-thermal case, using (9) for the remaining integrand, I obtain 

F L 2 < cos0' min [(v min /v ave )(\-l)/(\-2)] 1 -\ (15) 

For the L2 thermal case, I use (10) instead of (9), noting that the first two terms in (10) 
approximately cancel, giving 

F L2 ~ cos(j)' min {v min /v ave )Exp[-v 2 min /vl ve ]. (16) 

Now consider St = 61/ (\v\\\), the time between reflections, where SI is the length scale 
of the fluctuation 5B and (\v\\ |) is the average magnitude of the reflected particles' velocity 
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parallel to B T . For LI, (|v|||) = v\\ + /2 = (v ave /2)cos^ min ~ (v ave /2)(6B / B T ) X I 2 where 
the latter similarity follows for SB/Bt < 1 in (7). Thus 

5t Ll = (25l/v ave )(B T /5B) 1 / 2 . (17) 

For L2 ~ v\\ + /2 ~ (y A /2)cos 2 c/w ~ (V W /2)(5B/B T ), so that 

<ft L2 = (28l/V w )(B T /SB). (18) 

Consider now N(U/U r ) appearing in (8). For LI, U/U r ~ 1, as the energy of reflected 
particles is of order the average energy of all particles, but we must determine N. For LI, 
the energy gain is stochastic (Fermi 1949; Spitzer 1962; Eilek & Hughes 1991; Larosa et al 
1996) as the particles incur random walks through momentum space and 

N L1 = U 2 /(5U+} 2 , (19) 

where (5U + ), is the average energy gain by a particle from a head-on reflection. For L2, 
there are mainly head-on reflections, so that U r / (SU + ) < N L2 < U 2 / {5U + ) 2 . Since a lower 
bound on r r suffices, I employ Nl2 > U r / (8U+). This means that For L2, 

[NU/U r ] L2 >U/(5U+). (20) 

I now need (8U+) for both LI and L2. The (SU+) is determined by energy and 
momentum conservation before and after a mirroring. This gives {SU+) ~ 2mVA(\v\\ |). 
For LI, using the value of {\v\\\) calculated above then gives (5U+) ~ mVAV ave (5B / B) 1 / 2 , 
and thus 

[NU/Urhi = N L1 = {v 2 av jAV 2 ){B T /5B), (21) 
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while for L2, using the appropriate (|i>|||) calculated above (18) gives 

[NU/U r ] L2 >U/(6U+) = {v 2 ave /2Vl){B T /8B). (22) 

Collecting the calculations of N(U/U r ), St, and F for LI in (8) gives 

tlx = (St L1 /F L1 )[NU/U r ] L1 = (5l/AV A )(v ave /V A )(B T /5Bf = (l T /4V A )(v ave /V A )(5l/l T ) 1/2 i 

(23) 

where the last equality follows from assuming a Kraichnan (1965) spectrum (8B/B) = 
(51 /It) 1 / 4 relating the magnetic to scale fluctuations. Eq. (23) describes "stochastic 
Fermi" energization (Miller 1985; Melrose 1986). Similarly for L2, using the appropriate 
above results for N(U/U r ), St, and F in (8), I obtain 

t L 2 = {St L2 /F L2 )[NU/U r ] L2 >{5l/V A )(v av jV A )\B T /5Bf/ 2 Exp[Vl/vl ve ] 

= (lT/V A )(v ave /V A f(5l/l T f/ 8 Exp[V 2 /v 2 ave } : (24) 
for the thermal case, while for the non-thermal case with v m i n > vo, and A > 3 

tl2 > {lT/V A ){v ave /V A f- x {5l/l T f/\ (25) 

where the Kraichnan (1965) relation has again been used. 

We see that each of the energy doubling times (23)-(25) depend only the particle 
average speeds v ave and not on their particle mass. But if electrons and protons have the 
same v ave , the protons have (m p /m e ) more energy. Thus each of (23)-(25) shows that 
electrons take (m e /m p ) longer to drain the same amount of energy. When electrons and 
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protons do not have the same v ave one population could be in L2 and the other in LI and the 
comparison of energy doubling times becomes more subtle. This is because although in LI 
there are many more reflections possible than in L2 LI has both energy gaining and energy 
losing reflections (i.e. note the smaller shaded area and absence of symmetry in figure la 
compared to figure lb). Thus the energization is second order as expected for stochastic 
acceleration. For L2 however, while there are less reflections, they are mainly head-on 
(i.e. energy gaining). These two effects (less reflections but mainly energy gaining=Ll 
vs. more reflections but both energy gaining and energy losing =L2) compete and the 
j3 p regime for which protons vs. electrons dominate the drain then also depends on the 
particle distribution of that population in the L2 limit. Another complication comes if 
the populations have the same v ave but different distribution functions. Then one must 
compare (24) and (25). We will study some of these cases more specifically in the next 
section. 

5. Application to Accretion Flows 
5.1 Why Protons Are Likely Non-Thermal in ADAFs 

The discussion of section 4 is one approach to the mirroring or Fermi energization pro- 
cess. Others include stochastic magnetic pumping (e.g. Hall & Sturrock 1967) and transit 
time pumping (e.g. Stix 1962, described as the magnetic analogue of Landau damping). 
Achterberg (1981) showed that all small amplitude (5B « B T ) approaches in LI to mir- 
roring in a turbulent plasma can also be described by quasi-linear diffusion of particles in 
momentum space, from magnetosonic wave particle resonances at the Cherenkov resonance 
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(u> w — k\\vcos(j)) = 0. The relevant waves have frequencies u> w << particle gyrofrequencies 
(i.e. very long wavelengths compared to the gyro-radii) which is equivalent to the adiabatic 
approximation discussed in section 4. 

This resonance requires a minimum particle speed v m i n ~ Va and also a minimum 
sincf)' , as derived in section 4. The required minimum in sine/)' means that in order for 
particles to undergo repeated reflections and gain energy, their momentum must be rapidly 
isotropized on a time scale shorter than the time between reflections, which itself must 
be shorter than the largest eddy turnover time. The largest eddy turnover time is in 
turn shorter than the ADAF infall time, given by ti n ~ 1.8 x 10 -5 Mr 3 / 2 /a, where M 
is the mass in units of M , and r is the radius in Schwarzchild units. But for ADAFs, 
Coulomb isotropization is not fast enough: the time scale for momentum isotropization 
from Coloumb collisions is of order the time scale for thermalization and is given by (e.g. 
Spitzer 1960; Mahadevan & Quataert 1997) 

t pp = (27i) 1 / 2 (n p a T clnA)- 1 (m p /m e ) 2 (kT p /m p c 2 ) ~ 10" 2 a(/3 a /0.5) 3/2 MM, (26) 

where n p is the proton number density, ot is the Thomson cross section, Ink is the Coulomb 
logarithm, and M is the accretion rate is units of the Eddington value, 1.4 x 10 18 Mg/sec. 
Setting (26) equal to ti n shows that protons can only be Coulomb thermalized/isotropized 
well outside of the dominant energy emission location, i.e. for r > 100 (Mahadevan & 
Quataert 1997). (In fact, this feature is fundamental to enabling an ADAF solution.) Thus 
the isotropization requires an additional kind of wave-particle resonance. 

Unlike the mirroring waves, the required isotropizing waves have wavelengths of order 
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the particle gyro-radius. Some, or all these small wavelength (Whistler, Alfven or magne- 
tosonic) waves can be generated by the particles themselves and then they do not transfer 
energy to the particles. Some fraction may also be generated directly from the turbulence 
in which case they can transfer energy to the particles. This latter possibility is explored 
in (e.g. Quataert & Gruzinov 1998) as the primary means by which the turbulence dissi- 
pates into particles. The resonances occur when the wave frequency in the particle frame 
is an integer multiple of the particle gyrofrequency, that is ui — k\\vcos(f) — iVfi* = 0, where 
fi* = eBc/E p and E p is the total proton rest-f- kinetic energy. Quataert & Gruzinov (1998) 
show that Whistlers are not damped by protons. The short wavelength Alfven waves 
(u> = k\\VA < &> g = eB/(m,ic) < Q*) are the most relevant for isotropization and have 
the approximate resonance condition —eBc/E p — k\\vcos(f) = 0. The condition \cos(f)\ < 1 
then leads to the injection condition E p > (VA/v)miC 2 (Q g /uj). For u> ~ Q g , this leads to 
Vmin ~ Va for protons- similar to the requirement of the mirroring waves derived earlier. 

So both types of resonant waves-long wavelength mirroring waves and short wave- 
length Alfven waves (whether they accelerate or just isotropize)- have a proton minimum 
speed requirement of order Va- The inefficiency of Coulomb thermalization, and the need 
for wave particle resonances to dissipate the turbulence means that a significant non- 
thermal particle population should be produced. Since ADAFs are most commonly mod- 
eled with P a ~ 0.5, non-thermal protons will likely be a fraction ~0(1) of the population. 
The inefficiency of Coulomb collisions in ensuring a non-thermal population is fundamen- 
tal. Even when stochastic Fermi energization (the LI limit), can be shown to rigorously 
lead to a power-law distribution in the energized particles (e.g. Eilek & Hughes 1991) 
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efficient Coulomb scattering would thermalize the distribution. The fact that Coulomb 
collisions are inefficient, as shown above, precludes redistribution of energy over the full 
population of protons. 

Note that at least the isotropizing waves are also implicitly built into ADAFs because 
the standard models presume isotropic pressure and this would be impossible without wave- 
particle resonances. In fact, the plasma must be "collisional" in the sense of wave-particle 
interactions, even though it is "collisionless" with respect to Coulomb collisions. In short, 
non-thermal protons should be a generic prediction of ADAFs. This is consistent with 
observations (Mahadevan 1998), which can distinguish between thermal and non-thermal 
proton distributions in an ADAF framework (and so far are not too sensitive to the proton 
power law index.) In principle, similar arguments could be applied to electrons with more 
stringent resonance conditions. However Mahadevan & Quataert (1996) and Ghisselini 
et al. (1998) argue that synchrotron self-absorption can thermalize weakly relativistic 
and non-relativistic electrons (at least those not produced from pion decay) under ADAF 
conditions during an in-fall time. Thus I assume (25) applies to ADAF/HIT protons in 
the steady state, and (24) to electrons. 

5.3 Reflecting Waves, Scales of Dissipation, and when Mirroring Preferen- 
tially Energizes Protons vs. Electrons 

Magnetosonic waves are the dominant mirrorers in the low amplitude limit (SB << 
Bt), as Alfven waves are incompressible and compression is required for mirroring, though 
the relevant compression speed along the field lines is always ~ Va regardless of the 
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wave mode. Achterberg (1981) considered a magnetically dominated plasma at a single 
temperature, and focused only electrons. Here we are interested in a two-temperature 
plasma and are considering both electrons and ions. In general, both slow and fast waves 
may participate in the mirroring. 

Though magnetosonic waves dominate in the low amplitude limit, this is not necessar- 
ily true in the large amplitude limit (SB ~ B T ). Because the largest scales of turbulence 
in discs are the most compressive, the large amplitude limit is relevant when the the scale 
on which the mirroring can compete with the cascade of energy from larger to smaller 
turbulent scales is a large fraction of the outer turbulent scale It- In this case, the energy 
could be compressively drained into particles before it reaches smaller scales in the cascade 
where incompressible modes of dissipation dominate. Since large amplitude Alfven waves 
are compressible (Alfven & Falthammar 1963), even they could then contribute to the 
mirroring. Such Alfven waves could even steepen to form shocks and perhaps shock-Fermi 
acceleration would be relevant. This must be considered in future work, as will see that 
in fact the relevant mirroring scales can be large. 

I now proceed to estimate the scales on which the favored particles are energized 
and when protons vs. electrons are favored. For low luminosity sources, ADAFs require 
1 — fv < 0.01, implying an accretion efficiency < 1% of that for thin discs (e.g. Rees et 
al. 1982; Narayan & Yi 1995). The respective energization times then need to satisfy 

T p /r e <C = T p (l-fr,)/T e <10, (27) 

where the subscript p(e) indicates ions (electrons). Since the turbulence cascades from 
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large to small scales, I compare the scales of energy drain for protons, (51) p , and electrons, 
(51) e , for which (27) is satisfied. The larger of the two length scales then determines the 
dominant drain. Using (25) for ADAF protons and setting it equal to the eddy turnover 
time It/Va gives 

(5l) P /h ~ (V A /v p , ave )( 24 - 8X ^ 5 , (28) 

where v Pjave is the average proton speed. For electrons, setting ( times (24) equal to It/Va 
gives gives 

(5l) e /h ~ C 8/ HV A /v e , ave ) 2 ^Exp[-8V*/5vl ave }. (29) 
Then we can see that 

(5l) p /(5l) e = C 8 / 5 k 24 / 5 (3 4 p x / 5 Exp[8/(5k 2 (3 p )]. (30) 

This is > 1 for a range of parameters applicable to ADAFs (e.g. f3 p ~ k ~ 0(1), ( ~ 10). 
The same conclusion results when the particle distributions are either both thermal or non- 
thermal. Protons can be favored to the required extent when compressive modes dominate 
the turbulent dissipation. 

Let us determine the scale on which the protons are dissipated. From (28) we see 
that when V A ~ v ave , p , it is not hard to have (51) p /It ~ 0(1) > (recall A > 3). This 
means that the compressive modes may be very important and much of the energy in 
the turbulence may drain before approaching the incompressible scales where Quataert 
& Gruzinov (1998) is applicable. In general, the small amplitude limit may not be fully 
appropriate in describing energy dissipation in ADAFs. 
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Now consider a thermal plasma with T p = T e , which corresponds to radii outside 
the ADAF region (Narayan & Yi 1995) or to a thin, precursor disc. In this case, no 
matter which particles initially receive the energy, Coulomb collisions redistribute this 
energy between electrons and protons. However, whether protons vs. electrons receive the 
dissipated energy determines the heating rate. When j3 p ~ 0(1), the relevant limits of 
interest are (23) for electrons and (24) for protons. Using v e ^ ave = {m p /m e ) l / 2 v p ^ ave , the 
dissipation scale ratio becomes 

(Sl) p /(Sl) e ~ (m p /16m e )/3- 7 / 5 i?^[-8/(5/3 p )]. (31) 

This is < 1 for (3 P < 0.25 and > 1 for (3 P > 0.25. For (3 P » 1, both electrons and protons 
are in the limit of (23), for which (Sl) p /(Sl) e ~ m p /m e . For (3 P < m e /m p , both electrons 
and protons are in the limit of (24) for which (8l) p /(8l) e ~ (m p /m e ) 4 / 5 . In sum, electrons 
are favored only for the range m e /m p <(3 P < 0.25, while for j3 p outside this range protons 
are favored. (The conditions of low j3 p for which electrons are favored may be found in 
solar flares (e.g. Larosa et al. 1996) and some thin accretion disc coronae models (Field 
& Rogers 1993). 

5.4 Can Mirroring Help Form an ADAF? 

It is sometimes believed that purely a low enough accretion rate is enough to form 
an ADAF/HIT. However, unless the disc is already thick, the critical accretion rate below 
which electrons and protons do not couple by Coulomb collisions on an infall time as 
computed for a standard thin disc is far too low to be physically relevant. For a thin disc 
system to evolve into an ADAF, a mechanism is needed to form a thick disc first. This 
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may occur by thermal instability and mirroring may help. The condition (e.g. Pringle 
1981) for thermal instability is 



If the instability proceeds from within an optically thick disc, then we must compare 
blackbody emission to the heating. In the regime (3 P < 0.25 for the thermal disc, electrons 
are favored as shown in above, and (23) is applicable. Taking the inverse of (23) for 
electrons, multiplying by v\ ve oc T and differentiating gives dln[q t ]/dT = 1/2T. If the 
emission is blackbody, then 



and the instability is not favored. For (3 P > 0.25 protons are favored and using (24) gives 
dln[q t }/dT = T _1 (l//? p - 1/2T) and still, even for 0.2 < /3 p < 0.5 the thermal instability 
cannot ensue. 

However, it is more likely that the formation of an ADAF would proceed by thermal 
instability within the very surface layer of the thin disc, and successive layers would even- 
tually evaporate from the surface to form the thick ADAF disc. The particle distribution 
in the very surface layer could be non-thermal. To see how mirroring might help, in the 
limit where the protons dominate the energy drain ((3 P > 0.2), I invoke (25) for protons, 
take its reciprocal, and multiply by v\ ve oc T p = T to obtain the quantity proportional to 
q t . Then dln[q t ]/dT = (A — 2)/2T. For A > 8 this can satisfy (32) when the emission is 
blackbody (33). For Bremsstrahlung dln[q~]/dT = 3/2T, and (32) can be satisfied when 
A > 5. 



dln[q t }/dT > dln[q~]/dT. 



(32) 



dln[q-]/dT = 4/T, 



(33) 
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6. Conclusions 



Dissipation of turbulence in presumed ADAF sources must preferentially energize 
protons by a factor 77 > 99% over electrons if ADAFs are to account for the observed low 
luminosities. Compressive magnetic mirroring can in principle favor protons to the extent 
required for a ratio of particle to magnetic pressure (3 P > 0.25. This is less stringent than 
the requirements of incompressible modes which demand (3 P » 1. However, it is not 
easy to estimate what fraction of energy is dissipated in compressive vs. incompressible 
modes. An estimate of the scale at which mirroring can compete with the transfer of 
energy down the turbulent cascade seems to indicate that for j3 p ~ 1, the scale can be 
quite a large fraction of the outer turbulent scale — which for ADAFs/HIT can be a large 
fraction of the disc size (Blackman 1998). This has two implications: 1) a significant 
fraction of the energy may be dissipated in compressive modes and 2) the small amplitude 
approach to dissipation may not be valid. In the small amplitude limit, the relevant waves 
involved in the energy transfer to particles are magnetosonic waves, as Alfven waves are 
not compressive and will not be damped by mirroring. However on the larger scales in 
the turbulent cascade, the large amplitude limit is relevant. Since large amplitude Alfven 
waves are compressive, they too may be involved in compressive modes of dissipation. 
This presents additional complications for future work. Magnetic reconnection may also 
be a complication, however reconnection itself generates turbulence, and possibly shock or 
direct acceleration processes which may also favor protons, but it is important to know on 
what scale the reconnection is occurring. 
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Regardless of the fraction of energy dissipated in compressive or incompressible modes, 
the fact ADAFs are "Coulomb collisionless" on the radial infall time scale seems to make 
a non-thermal proton population inevitable. The required dissipation of turbulence must 
proceed through wave particle interactions, all of which act on only a subset of the particles. 
Since ADAF models presume an isotropic pressure tensor, wave-particle resonances are 
implicitly assumed to play a role in ADAFs because Coulomb isotropization is necessarily 
too slow. The presence of a non-thermal proton population seems to indeed be indicated 
by observations of the Galactic centre (Mahadevan 1998) when modeled with an ADAF. A 
remaining fundamental problem which still needs more attention is the question of a faster 
than Coulomb coupling between particles (Begelman & Chiueh 1988) even if the protons 
could receive the dissipated energy. 
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FIGURE 1 CAPTION 



Particle speed diagram for magnetic mirroring. The magnetic compression is assumed 
to move at velocity ~ — Va along the field line and so an isotropically distributed population 
of velocities in the lab frame has a spherical distribution centered around v'^ = Va- The 
angle (f)' min and the speed v m i n bound the respective minima needed for a particle to reflect 
at the magnetic compression. The weaker the compression the larger these minima. The 
area inside the shaded region between the two circles represents the particle speed region 
which can be reflected. Approximate schematics of the regimes a) L2 (v ave ~ Va) and b) 
LI (v ave » Va) of the text are shown. L2 is relevant for (3 a ~ 1 ADAFs. Note that only 
the region LI strictly corresponds to stochastic energization since in this case the number 
of "catch-up" and head on reflections are about equal whereas region L2 has primarily 
"head on" (albeit less total) reflections. 
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